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Abstract: We consider a holographic superconductor with homogeneous impurities added. 
We start with the holographic Abelian-Higgs model for s-wave superconductivity, and turn 
on a coupling between the gauge field and a new massive gauge field that is introduced 
for impurities, whose effect is examined in the probe limit. We find that the condensation 
of the massive gauge field is induced in the superconducting phase. When the coupling is 
sufficiently large, the mass gap in the optical conductivity disappears. A resonance peak 
is found in the conductivity for the massive vector field. 
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1 Introduction 

The AdS/CFT duality is a brilliant concept that allows us to describe strongly coupled 
conformal field theories with their gravity duals [1-3]. Certainties of this correspondence 
have been confirmed in so many examples, and it has been applied to various theories with 
much success. It is tempting to use this powerful tool to study strongly coupled systems 
with physical importance. 

Recently, applications to condensed matter physics have been intensively investigated. 1 
For example, superconductivity is modeled in gravity duals. First, it was pointed out that 
the black holes in AdS can accommodate scalar hairs [11]: a scalar field can have nonzero 
condensate outside the horizon of the black holes if temperature is sufficiently low. Based 
on this idea, holographic model of superconductor was then proposed in [12], and many 
generalizations were considered. (For instance, see [13-16]. See also references in [4- 
10].) Since gravity duals describe strongly coupled systems, it is hoped that holographic 
superconductors would unveil the mechanism of high-Tc superconductors. 

In condensed matter systems, considering the effect of impurities is often important 
since their presence can drastically change the physical properties of the systems. They can 
be used to control the conductivity of semiconductors. In superconductors, systems can 
show gapless superconductivity in the presence of impurities. Apart from such practical 
applications, it also generates physically interesting phenomena like Kondo effect. 

In this paper, we study effect of impurities for strongly coupled system using the gravity 
dual. We consider an s-wave holographic superconductor where two vector fields and 

are introduced in the AdS space: is dual to the conserved current in the boundary 
identified as the weakly-gauged electromagnetic field in the context of superconductors, 

1 For reviews, see [4-10]. 
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while the other is massive and dual to impurities. More concretely, we assume that the 
impurities have another type of charge carriers like a hole, for instance, where the hole 
number is not necessarily conserved due to the capturing of conduction electrons. So we 
assume that the latter is dual to the massive gauge field B^. We postulate a gauge invariant 
coupling term to describe the interaction between the two species. 

The model we consider was mentioned in [17], where the massive vector field could 
be introduced as impurities, but it was integrated out assuming that the impurities were 
infinitely heavy. If the mass is finite, the field at the asymptotic AdS boundary decays not 
exponentially but by a power law. Therefore, we do not integrate out the massive field, and 
solve the model explicitly in order to focus on the dynamics in the presence of the massive 
vector field. Our result will show that the mass gap of the superconductor disappears due 
to the interaction. 

The rest of this paper is organized as follows. In section 2, we study the model. The 
action and setups for numerical computations are provided in section 2.1. Numerical results 
are shown in section 2.2. Some further possibilities of the model are discussed in section 2.3. 
Comments on normal phase are given in section 3. In section 4, we reexamine the strategy 
of [17], and discuss relations to ours. We conclude this paper with future perspectives in 
section 5. 

2 Impurity degrees of freedom by a massive vector field 

We consider a holographic superconductor where an extra massive vector field is introduced 
as impurities. 

2.1 The model 

We consider a model where a massive vector field is introduced into the minimal Abelian- 
Higgs model of the s-wave holographic superconductor [12]. The action is 2 

S = j d 4 xv^( - ~F^ V ~ - iA^\ 2 - M 2 & 

- - A G, y G^ - —B^ - -F^j , (2.1) 

where F = dA and G = dB. The scalar field <3? is charged only under An, and there is 
no direct coupling between $ and B^. Here m 2 is the mass of B^. It might be possible 
to consider to generate this mass by some Higgs mechanism, but here we would like to 
start from the Proca action for simplicity. The scalar mass is chosen to be M 2 = — 2 for 
convenience in analysis. 

The action (2.1) has an interaction term between and G^. These vector fields are 
dual to two currents with different dimensionality, and this interaction represents a coupling 
of the two currents [17]. One current is conserved fermion number, and A^ is identified 
as the gauge field of the weakly- gauged U(l) electromagnetic symmetry in the context 

2 We set the gauge coupling of as e = 1. 
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of superconductors. The other, B^, represents impurities. If we allow the difference in 
anomalous dimensions of the two current operators, the vector field dual to the impurities 
can be massive and associated with the Proca field. 

We work in a limit in which the matter fields do not give back-reactions on the back- 
ground metric. 3 This limit is called a probe limit. The gravity background we consider is 
the AdS-Schwarzschild black hole, 

1 / dr 2 \ 

ds 2 = ^(-f(z)dt 2 + — + dx 2 + dy 2 ), f(z) = l-z\ (2.2) 

where we use units in which the AdS radius is unity, and the location of the black hole 
horizon is at z = 1. Such parameter fixings are possible thanks to the isometry of the AdS 
space, and can be confirmed by examining the symmetry of the equations of motion. 
The equations of motion of (2.1) are 

V A F A/i + cV A G v - 2|$| 2 ^ + i(<S>*d^<S> - <9 M <T$) = 0, (2.3) 
V A G V + cV x F Xl " - m 2 B^ = 0, (2.4) 

1 = (0„ - L4 M ) (^g^(d u - iA u )<S>) - M 2 <K = 0. (2.5) 



The kinetic terms of and B^ can be separated. From (2.3) and (2.4), we obtain 

(1 - c 2 )V A F A ^ - 2|$| 2 ^ + cm 2 B^ + i(^*d^^> - 5^*$) = 0, (2.6) 
(1 - c 2 )V A G A ^ - m 2 B^ + 2c\<5>\ 2 A^ - ic(^*d^ - d''$*$) = 0. (2.7) 

We can realize an superconductor phase where the scalar operator dual to $ acquires 
nonzero condensate (O). The ansatz is A = At(z)dt,B = Bt(z)dt, and = <p{z). The 
equations of motion become 

k2 A f cm 2 B t 



~ (l-c 2 )z 2 f + z 2 f 



,. m 2 B t 2cd> 2 A t 
z 2 f (1 - c 2 )z 2 f 



+ h7 + iH = °- ( 2 - 10 ) 



/' 2 \ jj ( 2 A t 

J--z)^ + \—f + T 2 

The prime denotes the derivative with respect to z: A' t = d z A±. We find it convenient to 
define an effective mass of B^, 



rn 



2 



m 2 = 2- ( 2 - n ) 

1 — cr 

The system is not well-defined if c = 1. If we consider a small deformation from a 
system without the impurity coupling, then < c < 1 would be reasonable. 4 For the 



3 If we recover the gauge coupling e and rescale each field by a factor of 1/e, then there is a factor 1/e 2 
appearing in front of the right hand side of (2.1). Backreactions on the metric are suppressed in the limit 
e 2 — > co with the fields fixed. The result takes the same form as (2.1), where e = 1. 

4 If c is negative, solutions of B M flips the sign. 
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strong coupling regime, we can take c > 1. In this case, however, At is tachyonic, and 
there can be some instability as we will discuss in section 2.3. We consider the case that 
< c < 1 first. 

Data of the boundary theory are extracted from the asymptotic behaviors of the bulk 
fields near the boundary (z — > 0). Let us start from c = case, where decouples from 
Ap and (f>. The asymptotic solutions to (2.8) and (2.10) take the form, 

A t = p-pz + ---, (2.12) 
<p = faz + (p 2 z 2 + ■ ■ ■ , (2.13) 

where p, p, (i = 1,2) are integral constants. Coefficients of higher-order terms are given 
in terms of these four constants. Here p and p are the chemical potential and charge density 
of the U(l) gauge field, respectively, and <f>i corresponds to the source or condensate of the 
scalar field. We impose boundary conditions such that the source of the scalar field is zero. 
We may consider either (f>± = or fa = 0: Since M 2 = — 2 is below the unitarity bound 
[18], the conformal dimension given by A(A — 3) = M 2 has two normalizable solutions 
A = 1, 2. The boundary conditions fa = and fa = correspond to A = 1 and A = 2, 
respectively. 

If c = 0, eq. (2.9) reduces to 

(2.14) 

z 2 f 

The asymptotic solution to this equation takes the form 

B t = f3„z a ~ + ••• + (/?+ + b + log z)z a + + ■■■ , a± = -(l±y/l + 4m 2 ) , (2.15) 

where /3± are integral constants: /3_ and /3+ correspond to the source and the condensate, 
respectively. The logarithmic term may appear only if a+ — a_ £ Z. The same statement 
applies later without mentioning explicitly. Solutions exist when the mass is above the BF 
bound [19], m 2 > —1/4. If m 2 = 0, the U{1) symmetry of B^ is restored, and a± = 0, 1. 

Now we take into account the coupling in the equations of motion. The series expan- 
sions (2.12), (2.13) and (2.15) are modified so that the asymptotic solutions accommodate 
the interaction term. For simplicity, we focus on the case that a± defined below are inte- 
gers. We also impose that the source of Bt is absent: /3_ = 0. The asymptotic solutions 
then take the form 

A t =p- pz + --- + (a+z a + +...) log* + ... , (2.16) 
B t = p + z a + + ■ ■ ■ + (b + z a + + ■ ■ - )logz + ■ ■ ■ , (2.17) 
= ((>iz + faz 2 + • • • + {<p+z a++2 + ■ ■ -)\ogz + • • • , (2.18) 



where 



a± = ^(l ± \A + 4m 2 ) . (2.19) 
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Note that the presence of the possible logarithmic terms is due to the impurity coupling. 
The expansions contain five integral constants fj,, p, (pi, and (p2- We further impose that 
the source of the scalar is absent: (pi = {i = 1 or 2). We will later show the asymptotic 
solutions in the cace that m 2 = 2. 5 

For the asymptotic forms of the solution at the horizon, we require that the gauge 
fields are zero at the horizon, while the scalar field can have a finite value. The series 
expansions take the form 

A t = A[ l) (z-1)+ A {2) (z - l) 2 + • • • , (2.20) 
Bt = B[ l \z-l) + B[ 2) (z-lf + ... , (2.21) 
cP = 0(°) + 0(1) (z-l)+ ( z - If + . . . , (2.22) 

where A^ , B^ and </>(°) are the integral constants, and higher-order coefficients are de- 
termined in terms of them. 

We will numerically solve (2.8), (2.9), and (2.10) to compute the condensate. Imposing 
the boundary conditions, we are left with four undetermined integral constants at the 
boundary, while there are three at the horizon. Thus we obtain one-parameter family 
solutions. Each solution is characterized by the temperature T of the system. 

It should be noted that Bt can become nontrivial in the bulk due to the coupling term 
in (2.9), which is nonzero if the scalar condensate is nonzero. In this case, Bt is inevitable 
to become nontrivial in the bulk even though the boundary conditions are Bt\ z= o tZ= ± = 
0. Eventually, we should obtain spontaneous condensation of Bt coming along with the 
condensation of the scalar field. 

Once we obtain nonzero condensate of the scalar field, we shall examine the proper- 
ties of this ordered phase by computing the optical conductivity. The ansatz of electric 
perturbations is A x = A x (z)e~ lu)t and B x = B x (z)e~ lu)t . The equations of motion are 

„// /' ,/ fu 2 2(P 2 \ A cm 2 B x 

+ + [p - (T^f/J x + ~z^r = ' (2 - 23) 

+ j B * + {p-^f) B * + O^Wf ~ °- (2 - 24) 

We impose the ingoing boundary condition at the horizon. The asymptotic solutions at 
the horizon hence take the form 

A x = r*"/ 3 (If + A^iz - 1) + A x 2 \z - I) 2 + • • • ) , (2.25) 
B x = r iuj/3 (S(°) + §P (z-1) + B x 2) (z - l) 2 + • • • ) . (2.26) 

Higher-order coefficients are determined in terms of the integral constants A^ and B x °\ 
One of these constants can be fixed thanks to the scaling symmetry of (2.23) and (2.24) 
under (A x , B x ) — > (XA X , XB X ) with some A. For instance, we may fix ^4^° = 1- The 

5 We obtain (a_,a+) = (-1, 2) when m 2 = 2. 
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asymptotic solutions at the boundary in general take the form 

A x = A x ~h a ~ + • • • + 4°) + A x l h + ■ ■ ■ + (5o + • • • + ~a+z a + + ■ ■ • )logz + • • • , (2.27) 
B x = B x ~h a ~ +■■■ + B x +h a + + ■ ■ ■ + (b + • • • + b + z a + + . . . ) logz + ■ ■ ■ . (2.28) 

There are four integral constants at the boundary: A x , A x l \ Bx ^ and B x + \ Other 
coefficients can be given in terms of them. 

The currents dual to these fluctuations A x and B x are computed from the on-shall 
action in the bulk. The action (2.1) can be separated into the bulk and boundary terms, 
S = Sbuik + where the former becomes zero because of the equations of motion, and the 
latter is used in computing correlation functions under the AdS/CFT prescription. The 
boundary action is obtained as 

Sb = \j d 3 x [A X 8 Z (A X + cB x ) + B x d z (B x + cA x )] z _> , (2.29) 

To compute the on-shell boundary action, we substitute the asymptotic expansions (2.27) 
and (2.28) into the boundary action, and then carry out the holographic renormalization 
[20, 21]. As a result, we obtain a quadratic action for fluctuations A x and B x . The 
quadratic action could be computed with general m 2 . For simplicity, however, we denote 

(2) 

the resultant quadratic action just as here, and will show an explicit computation in 
the case that m 2 = 2 in the next section. 

The AdS/CFT correspondence tells us that the currents of A x and B x are computed 
from the quadratic boundary action, 

A SS B SS 

= TnT' = TIT- (2.30) 

SA { X 0) 5B X } 

We call J x as an electric current, and J x as a hole current. Using these, we can give the 
two electric conductivity oa and 7# associated with two types of charge carriers, 



A 



cr A (u) = 3 -^- 
& x 



, tbH = -f- 



(2.31) 



where the external electric field applied on the system is given by E x = —dtA x °\ We 
numerically solve the equations of motion using the shooting method, and compute the 
conductivity in the superconductor background. 

2.2 Results 

For simplicity, we consider the case that m 2 = 2. We choose the coupling c = 0.5. The 
boundary condition of the scalar field is either 4>2 = or <f>i = 0, corresponding to A = 1 
or A = 2, respectively. In the figures, the orange dotted lines correspond to the case that 
c = 0, while the blue real lines correspond to the case of nonzero c. 
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First, we compute the condensate for superconductor background. As announced, we 
show the asymptotic solution at the boundary when m 2 = 2, after /3_ = is imposed, 

A t = pt - pz + L<ft - c/3+ + 3 2 5^ ^ log z^j z 2 + --- , (2.32) 

H^-^H' 2+ C ^(i-T 2) ' 3+ ---- <2 - 33) 

<p = 4>iz + (f> 2 z — z H z*H . (2.34) 

2 D 

Figure 1 shows the results of scalar condensate. The critical temperature Tq does not 
change regardless of c because Bf = at T = Tq, where (ft = 0. Thus we normalize the 
temperature by Tq ~ 0.23^fp in the case that A = 1, and Tq ~ 0.12y / p in the case that 
A = 2. We find that the scalar condensate is almost the same as the c = case, although 
there can be seen slight deviations in low temperatures. 

Figure 2 shows the condensate of B t , given by (£>) = (1 — c 2 )/3 + + c/x0f/3. The 
absolute value of (£>) is plotted with a power of 1/3. 6 The source of Bt is absent, and 
the condensation is triggered by /i, but the condensation is present only if the scalar field 
has nonzero condensate. Figure 2 tells us that the qualitative behavior of (B) is similar 
to that of {O). However, {B) goes to zero linearly at T — > Tc- It looks diverging in low 
temperatures in the case that A = 1, and backreactions on the gravity background will be 
important in low temperatures. The condensate (B) becomes nonzero once c is nonzero, 
and the magnitude becomes larger as we increase c. 

We compute the conductivity in the superconductor background obtained. First, we 
prepare the explicit forms of ct(uj) and 7(w) in the case that m 2 = 2. The asymptotic 



3 The sign of (B) is flipped if the sign of c is negative. 
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I<Ba=i>I 1/3 



I<Ba=2>|' /3 




Figure 2. The condensate of B t as a function of temperature when m 2 = 2, c = 0.5. In both 
panels, the condensate goes to zero linearly at Tc'- (B) oc T — Tc- 



solutions of A x and B x at the boundary (2.27) and (2.28) are obtained as 
A x = - ^ + 4 0) + f 4 } " 2c log * 



+ U 0) (0 



+ 



2c 2 4°Vi + 2 Cj Bi- } (l - c 2 (l + 4010a)) 



3(1 - c 2 ) 



log z \z + 



B 



(-) 



+ 



2 i- c y x 

i(0)^2 , odH 



2cA!b u V! + 2i?i" ; (l - c 2 (l + 40i0 2 )) 



log z ]z + 



3(1 - c 2 ) 

By making use of them, the currents dual to A x and B x are computed as 



J A 



6S, 



(2) 



4D + cflH " 



5c 2 )0 2 



3(1 -c 2 ) y 



J" 



_ 55 b 



(2) 



(1 - c 2 )4+) + c4 0) ^ 



a; 2 (8 - 5c 2 )0| 



2 3(1 - c 2 ) 



Thus, two conductivity are given by 
g a (uj) -- 

7b M = 



w 4 0) 



4 '=0 



w 4 0) 



(i- c 2 )4 +) + c4 0) (y-^ ( 



5c 2 )0 2 



4 '=o 



(2.35) 



(2.36) 



(2.37) 
(2.38) 



(2.39) 
(2.40) 



Notice that the second term in 7b can be found only in the imaginary part. 

In figure 3, we plot the real part of the electric conductivity as a function of ui normal- 
ized by the scalar condensate. The conductivity is computed for T /Tc = 0.60, 0.30, 0.20 
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2 - 5 V<o 2 > 



Figure 3. The real part of the electric conductivity as a function of uj normalized by the condensate. 
The left panel is when A = 1, and the right panel is when A = 2. The blue lines correspond to 
c = 0.5. The orange lines go to zero exponentially in low temperatures. However, the blue lines 
follow power-law behaviors even in low temperatures. 




2.0 




2.5<0|> 



2 - 5 V<o 2 > 



Figure 4. The imaginary part of the electric conductivity as a function of uj normalized by the 
condensate. The left panel is when A = 1, and the right panel is when A = 2. 



when A = 1, while T/Tq = 0.70, 0.35, 0.20 when A = 2. The rightmost lines correspond 
to the lowest temperatures. 

Let us first recall the interpretation in the case that c = (orange lines). In the 
superconducting phase, there is a particular feature that a gap turns up when T/Tq is 
lowered. There is no gap until some T = T g as Re(<r)| iJ _>o is wen above zero. However, 
when T/Tq is small so that Re(o")| w _j.o is close to zero, the gap quickly shows up when 
T <T g . This suggests that no quasiparticles are excited up to some uj. 

If there are the coupling to the impurities, however, agitation appears in low frequencies 
as observed in figure 3. The exponential growth of the conductivity in the case of c = is 
now replaced with a power-law behavior. The curve in the small frequency region can be 
fitted with a power function, 



Re(cr) oc 



(2.41) 



where Ai = (Oi) and A2 = \fIP2). Thus, this phase is gapless. 

Figure 4 is the imaginary part the conductivity as a function of uj normalized by the 
scalar condensate. Im(cr) is diverging as u — > 0. This behavior corresponds to the presence 
of the delta function at uj = of the real part of the conductivity. The plot lines slightly 
shift toward small-o; direction. 



- 9 - 




Figure 5. Comparisons of the real part of the conductivity for c = (orange dotted lines) and 
c = 0.2, 0.3, 0.4, 0.5 (blue real lines) computed when T/T c = 0.20. The left panel is when A = 1, 
and the right panel is when A = 2. Exponential and power-law behaviors are interpolated by 
changing c. 



It is expected that the effects of the coupling turn up gradually as the coupling is 
increased from c = 0. To discuss this, we compare cases of different c. We compute the 
conductivity for c = 0.2, 0.3, 0.4, 0.5, by fixing fh 2 for simplicity. We have to be careful 
that this gives the results corresponding to different values of m 2 in the action (2.1). For 
instance, if m 2 = 2 is fixed, the computations are for m 2 = 1.92, 1.82, 1.68, 1.5 when 
c = 0.2, 0.3, 0.4, 0.5, respectively. 7 However, we see that comparing with fixed fh 2 can 
give qualitatively reasonable understandings for the behaviors of the conductivity under 
change of c. In particular, we can discuss the presence or absence of the mass gap. 

Comparisons of the real part of the conductivity with respect to different c are given 
in figure 5. The graphs are computed when T/Tq = 0.20. We see that, as c is increased, 
the conductivity gradually changes from exponential to power-law behaviors in the small 
frequency region. It looks that the transition is continuous, and there would not be a clear 
phase transition from gapped to ungapped phases under the increase of c. 

In figure 6, we would like to present a schematic description of the phase diagram of 
our model under the change of c. The strength of the coupling c might be interpreted as 
the density of impurities. We expect that the regions of the gapped and ungapped phases 
are smoothly connected under the change of the density. The boarder would not be given 
by a clear phase transition. When c = 0, the gap starts to show up at about some T g . 
As c increases, the temperature for the appearance of the gap looks to be lowered. There 
may be some critical c = c* above which there is no gapped phase even at T = 0, or c* 
may be located at c — > 1. It will be interesting to look for such a quantum critical point. 
However, to discuss the zero temperature limit, we will need to carry out precise analysis 
by including back reactions on the gravity backgrounds [16]. 

We may consider to vary the mass fh 2 . When this is heavier, we find that the effects 
on the conductivity is smaller. Technically, this would be because a+ is larger and 
decays faster at the boundary. Hence the effects on A x from B x would be smaller at the 
boundary. If the mass is heavier, the massive vector field becomes more non-dynamical, 

7 If we fixed m 2 = 1.5, the effective masses would be m 2 = 1.56, 1.65, 1.79, 2 when c = 0.2, 0.3, 0.4, 0.5, 
respectively. We would do computations by employing different asymptotic expansions with respect to 
different values of fh 2 . 
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(O) = 0, (B) = 



(O) £ 0, (B) + 



gapped 

Re(cr) ~ e u 



ungapped 

Re(cr) ~ lu 2 
interpolating 



Figure 6. A schematic description of the phase diagram expected. Here (B) = if c = 0, while 
(B) ^ once c ^ and T < Tc- Gapped and ungapped phases would be interpolated smoothly. 
There may be some critical c = c* above which there is no mass gap. Figure 5 corresponds to 
looking at a horizontal slice of this diagram at T /Tc = 0.20. 



and thus the effects might be harder to appear. 

Let us come back to the case that m 2 = 2 and c = 0.5. The results of 7(0*) are 
given in figures 7 and 8, corresponding to the real and imaginary parts, respectively. The 
conductivity is zero at T = Tq since B x = in the normal phase, although we do not plot 
T = Tc case in the figures. As the temperature is decreased from T = Tc, the conductivity 
starts to grow. There cannot be seen a gap in the small frequency region of the real part 
as in the case of the electric conductivity. Hence, Re(7) can be fitted with a power-law 
function, 

ReH WfV, (2-42) 



where Ai = (Oi) and A2 = y/ (O2), and Ro = lim aJ _ s .o Re(7)/Aj is an intercept. In figure 7, 
a peak can be found in Re(7) at w\ ~ 1.2 when A = 1, and 102 ~ 1 when A = 2. The peak 
becomes more apparent as the temperature is lowered, and it is expected that Re(7)/Aj 
converges to a finite height curve in the low temperature limit. The nonzero slope in the 
large-w region of Im(7) is due to the direct contribution of the electric field. 

Similar to the case of the electric conductivity, we compare the dependence of 7 on c. 
The results for the real part when T /Tc = 0.20 are given in figure 9. Since is decoupled 
from when c = 0, 7 is not induced by the electric perturbations in this limit. From 
this figure, it can be seen that the position of the peak slightly moves to larger oj as c is 
increased. The magnitude of the peak looks to become larger as c is increased. However, 
since Re(7) is proportional to (1 — c 2 ), the magnitude of 7 would be suppressed if c is 
sufficiently large. It is expected in the heavy mass limit of m 2 that would be hard to 
be excited, and this would mean suppression of 7 in the c — > 1 limit. 
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Re(y) 




Re( 7 ) 



2.5 <Oi> 




2-5 V<0 2 > 



Figure 7. The real part of 7 as a function of us normalized by the condensate. The left panel 
is when A = 1, and the right panel is when A = 2. Larger magnitude lines correspond to lower 
temperatures. 



Im(y) 




Im(y) 

V<o 2 > 
0^ 



.5<Oi> -5 



jxi — iifi 11 -n 



2> 



Figure 8. The imaginary parts of 7 as a function of ui normalized by the condensate. The left 
panel is when A = 1, and the right panel is when A = 2. Larger magnitude lines correspond to 
lower temperatures. 



Re(y) 




Re(y) 

V<o 2 > 

0.25 



2.0<O!> 




Figure 9. Comparisons of the real part of 7 for c = (orange dotted lines) and c = 0.2, 0.3, 0.4, 0.5 
(blue real lines) computed when T/Tc = 0.20. The left panel is when A = 1, and the right panel 
is when A = 2. Lines with larger magnitude correspond to larger c. 



2.3 The possibility of negative effective mass 

In the previous section, we considered the case that < c < 1 and m 2 = 2. However, 
since the BF bound for the massive vector field is m 2 = —1/4, it can be possible that the 
effective mass of B^, m 2 = m /(l — c 2 ), is negative. One possibility is that the system is 
in strong coupling: c > 1 with m 2 > 0. The other possibility is that the mass squared is 
negative m 2 < but < c < 1. In the latter case, we obtain \m 2 \ < \m 2 \. Therefore, m 2 
satisfies the BF and unitarity bounds, —1/4 < m 2 < 3/4, if in 2 does. We discuss these 
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Figure 10. The conductivity of A = 1 case when m 2 = —1/4. The blue real lines correspond to 
this negative effective mass case with c = 0.2, while the orange dotted lines are when c = 0. The 
real part of the conductivity increases near w = 0, but there is a delta function at to = 0. 

options. 

First, we consider the possibility if the system is in strong coupling. In this case, we 
achieve 1 — c 2 < 0. Hence, the effective mass squared of At in (2.8) becomes negative. It 
is therefore concerned that the effective mass function of At may violate the BF bound at 
some point where (j) is large, although m 2 for Bt is above the BF bound. In fact, numerical 
examinations in the probe limit tell us that At can be badly tachyonic. Therefore, the 
assumption of strong coupling seems unreasonable. 

Second, we discuss the case that m 2 < but < c < 1. The scalar condensate is almost 
the same as that found in figure 1. Moreover, (B) goes to zero linearly at T — > Tc as well, 
although the slope can be different from the previous case. The conductivity, however, 
behaves strangely in this case. Figure 10 shows the electric conductivity when m 2 = — 1/4 
and c = 0.2. Curiously, the real part grows up as u decreases. This may look similar to 
the presence of the Drude peak. However, as is suggested by the imaginary part, there is 
a delta function at u = of the real part. It is certain that this rise in small frequency is 
due to the coupling to B x . In this case of c = 0.2, the mass gap looks disappearing owing 
to the uplift of the conductivity. It will be interesting if relevance of these behaviors to 
real-world condensed matter systems are understood. 

3 Without scalar field: RN-AdS black hole background 

We have considered the addition of the massive vector field in superconductor phase. How- 
ever, it will be interesting if we can observe some effects in normal phase. We try a case 
that the gravity background is the RN-AdS black hole, where the backreaction of the U{1) 
gauge field is included. Since the massive vector field is not charged, we do not turn 
on Bt, and hence the gravity background is not affected by B^. Here we discuss that 
spontaneous induction of B^ do not occur. This suggests that simply using the RN-AdS 
black hole might be too naive. It would be interesting to consider more general black hole 
backgrounds where B^ can back-react on the gravity backgrounds. 

Let us consider a model that contains a U(l) gauge field and a massive vector field 
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but no scalar field. The action is 



By denning the effective mass as m 2 = m 2 /(l — c 2 ), we obtain the equations of motion of 
the matter fields as 

V x F Xfl + cm 2 B^ = 0, (3.2) 
V x G Xfl - m 2 B» = 0. (3.3) 

These should be accompanied with the Einstein equation. 

We try the the Reissner-Nordstrom AdS (RN-AdS) black hole as the gravity back- 
ground, 

1 ( dz 2 \ 

ds 2 = ^ (-f(z)dt 2 + — + dx 2 + dy 2 \ , (3.4) 

/(*) = 1 + (m/7) V " (1 + W7) V, (3.5) 
A t = M (l - z), (3.6) 

where [i is the chemical potential, and /x/7 gives the electric charge density of the black 
hole. The location of the outer horizon is at z = 1. The Einstein equation can be solved 
with B^ = 0. Note that B^ is decoupled in (3.3). Without source term, B^ = can be the 
solution. 

Let us consider the optical conductivity of A x in this background. The ansatz of 
electric perturbation is A x = A x (z)e~ lujt , B x = B x (z)e~ tujt , and gt x = gtx{z)e~ lujt . The 
linearized equations of motion for the gauge fields are 

,// /' ,/ „ z 2 A' t ( , 2 \ cm 2 

A" x + J -A' x + —A x + — 1 (g' tx + -g tx + -—B x = 0, (3.7) 



/ * \f 2 z 2 f 
The linearized Einstein equation gives a first-order equation of gt 



2 4j4' 

<4 + ~9tx + ^{A X + cB x ) = 0. (3.9) 

z 7^ 



We can eliminate gt x from (3.7) by using (3.9). We obtain 



„„ /' (u? Az 2 (A') 2 \ A ( m 2 Az 2 (A') 2 x 

At the horizon, we impose the ingoing boundary condition such that A x and B x are pro- 
portional to f- tu) /(3-v 2 'hi 2 ) _ Unfortunately, (3.8) is decoupled from A x , and the solution is 
B T = 0. Thus, there is no contribution of B~. 
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4 Comments on adding impurities by bulk source terms 



In previous sections, we considered a model where there is a dynamical massive vector field. 
On the other hand, an idea on introducing impurities in holography non-dynamically was 
proposed in [17]. Since our model is related to theirs, we would like to discuss [17] in this 
section. 

The proposal is that one can introduce impurities by adding a bulk source term, 
J d^xA^J^, where the bulk U(l) source would correspond to the impurities. La- 
grangians with this bulk source term is treated as starting places for analyzing impure 
systems. This coupling is considered in the s-wave holographic superconductor. The ac- 
tion is 

S = f d 4 x^=g {- l -F^ u F^ - \d^ - iA^\ 2 - M 2 <5> 2 + A^jA , (4.1) 

where M 2 = —2. The probe limit is assumed, and the background is the AdS-Schwarzschild 
black hole (2.2). Here we reexamine this model. In particular, we shall understand how 
the effects of appear more precisely than the results in [17]. 

Taking (4.1) as the starting place may be regarded as an economical simplification in 
analyzing impurity systems. In giving (4.1), it is assumed that the impurities which would 
have been introduced in underlying theories are somehow replaced with J^. We may 
consider many possible Lagrangians which might reduce to the form of (4.1). For instance, 
the massive vector model that we focus on in this paper can be found in [17], where the 
massive vector field is interpreted as impurities. However, it might be basically hard to 
derive by integrating out the impurity fields explicitly. Therefore, one would start from 
an effective theory where the bulk source term is introduced as in (4.1). However, the form 
of should be arbitrary, and therefore needs to be chosen by hand. Some reasonable 
may be assumed. 

An example considered in [17] is to add impurity density by turning on J*, while other 
components of are zero. Using the ansatz A = At(z)dt and <J> = (j)(z), we obtain the 
equations of motion for scalar condensate, 



f+(^--)^+(4i + ^)<l> = 0. (4.3) 



2cj) 2 A t _ J 1 

v/ z)* ' \z 2 f ' f 2 

Since J* needs to be chosen by hand, J* = is adopted. This is devised such that J* 
decays fast enough at the boundary. However, how to choose J 1 is irrelevant to dynamics 
of the system. In particular, in this example J* 7^ at T = Tq- This is different from what 
we observed in section 2, where is zero at T = Tq, and nontrivial when T < Tc- 

Since J x is not present, the equation of motion for the electric perturbation A x = 
A x (z)e~ tu)t is as usual, 

K + f jA! x + A x = 0. (4.4) 
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c 


-1 


-0.5 


-0.1 





0.1 


0.5 


1 


Tc/T*\a=i 


1.058 


1.028 


1.005 


1 


0.995 


0.975 


0.951 


Tc/T* A=2 


1.020 


1.010 


1.002 


1 


0.998 


0.991 


0.982 



Table 1. Difference of the criticai temperatures from the c = case. Here T» is the critical 
temperature when c = 0. 



K + T4 " ^ + 4r = 0, (4.5) 



It is expected that the impurity density J* will give some effects on (ft through At, and the 
conductivity will receive contributions of J* through (p. Results obtained in [17] deduce 
that, compared at some given temperature, the real part of the conductivity increases in 
the presence of J*. 8 This would sound similar to our results. Hence, we would like to find 
out the grounds that lead us to this observation. 

The case of nontrivial J y is also considered in [17]. The equations of motion for the 
superconductor background need non-zero A y - component: 

f 2^ A y jy 

f y z 2 f z*f 

together with equations for At and (p. Here J* = 0. However, the conductivity equation 
takes the same form as (4.4). The upshot is that the results of the conductivity are 
qualitatively the same as the case of the impurity density J*. For this reason, we discuss 
only the case of the impurity density. 

We compute the scalar condensate and critical temperatures for superconductivity for 
various c. We use the same impurity density as in [17]: ,P = cz & . Note that the relative 
signs between c and fj, are important, and \x is positive in our convention. We find that 
the critical temperature Tq shifts depending on the magnitude of c. Results are given 
in table 1. For this reason, it is convenient to normalize temperatures by the critical 
temperature when c = 0, which we define as T*. We see that Tq decreases if c > 0, while 
increases if c < 0. The shifts are larger when A = 1 than when A = 2. 

Figure 11 compares the scalar condensate when c = 0, ±0.5. Effects of J* are recog- 
nized only near T#, and disappear when T/T* is small. It would be reasonable to think 
that the addition of impurities change the critical temperatures for the superconductivity. 
The reason why such shifts of the critical temperature occur seems to be because J* 7^ 
at T = Tq, where it looks that the configuration of At at (f> — > is affected by J*. This 
may be interpreted such that external operations are applied on the system through J*. 

It is expected from figure 11 that the effects of J* on the conductivity are significant 
only when T/T* is close to 1. Hence, we would like to look into this region. We compute 
the conductivity when T/T* = 0.95, and the results are shown in figures 12 and 13. We 
find that, in low frequencies, the real part of the conductivity increases if c > 0. This 
agrees with the results of [17]. However, we also find that the conductivity decreases if 
c < 0. The imaginary part decreases if c > 0, and increases if c < 0. These shifts are due 
to the difference of 6. 



8 A positive c was used in their computation. To be precise, the signs of c and chemical potential /i are 
the same. Conventions are such that the chemical potential /i is also positive. 
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Figure 11. The scalar condensate for A = 1 (left panel) and A = 2 (right panel). The blue real 
lines are when c = 0.5, the orange dotted lines are when c = 0, and the purple dashed lines are when 
c = —0.5. The temperature and condensate are measured in units of T*. The critical temperatures 
arc different for different values of c. 

Re(cr) 
1.0 

0.8 

0.6 

0.4 

0.2 



10 



15 



20 t 



Im(o-) 


3.0 


1 




1 

1 


2.5 


1 
1 
1 


2.0 


1 
1 

1 


1.5 


1 
1 
1 


1.0 


t 
I 
\ 


I 


0.5 


\ 

\ *,\ 






0.0 





10 



15 



20 t 



Figure 12. The conductivity when A = 1 and T/T* = 0.95. The blue real lines are when c = 0.5, 
the orange dotted lines are when c = 0, and the purple dashed lines are when c = —0.5. 
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Figure 13. The conductivity when A = 2 and T/T* = 0.95. The blue real lines are when c = 0.5, 
the orange dotted lines are when c = 0, and the purple dashed lines are when c = —0.5. 



We thus find that it is simply the shift of Tc that is relevant for the shifts of the 
conductivity in this model. This can be understood as follows: Owing to the shift of Tc, the 
magnitude of (ft at some T/T* changes, and this results in the shift of the conductivity. The 
effects of this shift of Tq disappear when T/T* is lowered, as differences in the condensate 
disappear. The mass gap in the conductivity appears in low temperatures. This is the 
same as the case that there is no impurity. 

This result suggests that there is no dynamical effects of impurities in low temperature 
in this model. Indeed, the impurities in (4.1) is non-dynamical. The effects on the conduc- 
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tivity are present only near T*, and this is due to the shift of Tq- On the other hand, the 
action (2.1) includes the dynamical massive vector field, and there are nontrivial results in 
low temperatures. It should be noted that Tq does not change in our model. 

5 Conclusion and discussions 

In this paper, we considered a model of an s-wave holographic superconductor where a 
massive vector field for representing impurities was introduced. The model had an coupling 
of the gauge field and the massive vector field B^, and was analyzed in the probe 
limit. We found that the massive gauge field was excited in the superconducting phase. A 
schematic phase diagram expected from our examination was given in figure 6. We also 
computed the optical conductivity for electromagnetic perturbation. When the coupling 
was sufficiently large, the mass gap in the conductivity disappeared. A resonance peak was 
found in the real part of the conductivity 75. We also considered the case of normal phase 
in the RN-AdS black hole background, and found there was no effect on the conductivity. 
We discussed the model of [17], where impurities were proposed to be introduced non- 
dynamically by adding a bulk source term. We have shown that, in that model, the critical 
temperature for the superconductor phase were shifted due to the presence of the bulk 
source. 

There are many things to be considered in the future. Several questions to be asked 
can be as follows. Perhaps the most straightforward extension of our work is to include 
backreactions of the matter fields to gravity backgrounds, as considered in [15]. We may also 
be able to consider non-relativistic deformations like the Lifshitz geometry [22, 23], which 
would be interesting as application to real-world condensed matter physics. We may also 
apply magnetic field to our model, and see the fate of the Meissner effect. The interaction 
we considered, cF^ w G^ u , is added in a bottom up way. Since this interaction is crucial for 
the results in this paper, it would be interesting to derive this interaction from a top-down 
approach in string theory. This interaction might appear in non-linear actions of D-branes. 
(Some related discussions can be found in [24].) Other interesting situations are to allow 
spatially inhomogeneous configuration as recently studied. In particular, realizing lattice 
structure [25, 26] is an interesting approach in modeling real- world materials. Another 
interesting case is to realize Josephson junctions [27]. We hope to come back to report on 
these topics. 
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